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Transient  Refractive  Index  in  a  XeF  Laser 

Distinguished  Professor  Allen  E.  Fuhs 
Department  of  Physics  and  Chemistry 
and  Department  of  Aeronautics 
Naval  Postgraduate  School,  Monterey,  California  93940 

CAPTAIN  James  Etchechury,  USA,  Graduate  Student 
LT  Lonnie  Cole,  USN,  Graduate  Student 
Naval  Postgraduate  School,  Monterey,  California  93940 

Gregory  A.  Blaisdell,  Graduate  Student 
California  Institute  of  Technology,  Pasadena,  California  91125 

Transient  refractive  index  was  calculated  for  an  e-beam 
pumped  XeF  laser  having  initial  concentrations  Ne:Xe:NF  :: 
95.5:4.3:0.2.   The  calculation  used  Ne  in  excited  states  in- 
cluding the  transition  array  3s  ->  4p  and  all  other  constituents 
in  ground  state.   Lines  in  3s  ■*■  4p  array  may  be  resonant  with 
XeF  radiation.   To  obtain  the  transient  populations  for  Ne,  a 
modification  was  made  to  a  Naval  Research  Laboratory  computer 
code  which  typically  accounts  for  180-200  reactions  involving 
50-60  different  species.   The  modified  program  included  popu- 
lations in  12  electronic  configurations  of  Ne.   Electron  beam 
current  was  changed  10%  at  constant  voltage  to  simulate  the 
influence  of  a  hibachi  shadow.   For  a  nonresonant  laser  wave- 
length, an  optical  pathlength  in  the  laser,  L,  less  than  35 
meters  will  have  distortion  6/X  less  than  0.1;  however,  for  a 
laser  wavelength  in  resonance  with  Neon,  L  of  only  0.11  meter 
will  cause  6/X  =  0.1.   The  difference  in  optical  pathlength 
is  6. 


INTRODUCTION 

Excimer  lasers,  which  provide  radiation  in  the  visible  and 

ultraviolet  regions  of  the  spectrum,  are  scalable  to  high 

1-3 
energy.     Because  high  energy  per  pulse  is  available,  excimer 

lasers  have  many  potential  applications.  Several  applications 

require  excellent  beam  quality. 

As  with  many  gas  lasers,  beam  quality  depends  strongly  on 

4 
the  refractive  index  field  within  the  lasing  gas.   The  pro- 
cesses of  pumping  and  power  extraction  may  generate  spatial 
gradients  of  the  index  of  refraction.   Further,  high  energy 
lasers  employ  flowing  gases  as  a  means  of  eliminating  heat; 
the  volumetric  release  of  heat  is  a  consequence  of  laser  in- 
efficiency.  Excimer  lasers  circulate  the  gases  which 
constitute  the  laser  medium. 

Excimer  lasers  are  inherently  a  pulsed  device.  When  the 
energy  is  deposited  by  the  e-beam,  the  static  temperature  of 
the  gas  may  increase  by  30%,  driving  strong  acoustic  waves  and 
even  shock  waves  with  Mach  number  the  order  of  1.15.   The 
impact  of  the  waves  on  the  index  of  refraction  must  be  elimi- 
nated between  pulses.   Acoustic  waves  are  phenomena  which 

7  8 
excimer  lasers  share  with  other  pulsed  electrical  lasers.  ' 

For  lasers  operating  with  vibrational  transitions  in  the 

infrared,  e.g.,  C02  and  DF,  the  index  of  refraction  can  be 

modeled  adequately  in  terms  of  the  Gladstone-Dale  constant  and 


the  variation  in  mass  density.   For  an  excimer  laser,  the  re- 
fractive index  does  not  depend  solely  on  mass  density.   As  the 
discussion  in  Section  I  demonstrates,  the  index  of  refraction 
depends  on  both  the  mass  fraction  of  chemical  species  and  the 
populations  for  each  of  the  chemical  species  in  addition  to 
mass  density.   Hence,  to  understand  variations  in  index  of 
refraction,  the  variations  in  chemical  species  and  populations 
must  be  understood.   Section  I  also  discusses  the  methods  to 
calculate  the  index  of  refraction. 

Variations  in  the  index  of  refraction  cause  distortion  of 

9  10 
phase  at  the  output  aperture  of  the  laser.  '    A  nonuniform 

distribution  of  phase  causes  a  severe  degradation  in  beam 

quality.   One  approach  to  calculating  the  index  of  refraction 

and  its  impact  on  beam  quality  is  to  use  the  Ladenburg  formula 

which  is  based  on  the  Kramers'  quantum-theoretical  dispersion 

formula.    The  Ladenburg  formula  requires  extensive  knowledge 

of  oscillator  strengths  and  photoionization  cross  sections  for 

molecules  in  the  ground  state  as  well  as  in  excited  states. 

Sources  of  data  are  discussed  in  Section  II.   Additional 

knowledge  concerning  the  rise  and  fall  of  chemical  species  and 

the  populations  also  is  needed.   A  computer  program  yielding 

the  information  on  species  concentrations  and  populations 

is  discussed  in  Section  II  and  III. 


Results  of  calculations  of  the  index  of  refraction  are  pre- 
sented in  Section  IV.   Combining  the  information  from  Sections 
III  and  IV,  an  estimate  of  the  influence  of  a  spatial  varia- 
tion in  input  power  density  on  beam  quality  is  made  in 
Section  V. 

I.   INDEX  OF  REFRACTION 
A.   Calculation  of  index  of  refraction 

The  development  in  this  section  follows  closely  that  of 

12 
Fano  and  Cooper.    The  dielectric  constant,  e;  complex  index 

of  refraction,  n;  and  electric  susceptibility,  x  »  are  related 

by 

£  =  n2  -  (n  +  iK)2  =  1  +  4irx  (1) 

The  real  part  of  n  is  of  interest  and  is  given  by 

1/2 
n  =  Re[l  +  4ttx   ]  =  Re[l  +  2ttx   ]  (2) 

Fano  and  Cooper  use  the  susceptibility  in  the  form 

S       OJ       -    00       -    iY    OJ 

s  's 

where  e  is  the  charge  of  an  electron,  m  is  mass  of  an  electron, 

N  is  the  number  density  of  molecules,  f  is  the  absorption 

oscillator  strength,  oo  is  the  frequency,  and  y     is  the  damp- 

s  s 

ing.   Further,  Fano  and  Cooper  point  out  that  when  photo- 
ionization  processes  occur,  the  form  of  Eq.  (3)  should  be 


discrete     °°  ,  , 

s       s       I 
In  Eq.  (4),  I  is  the  lowest  ionization  threshold.   Fano  and 
Cooper  discuss  the  influence  of  autoionization. 

A  relation  exists  between  the  photoionization  cross  section 

and  the  spectral  density  of  oscillator  strengths  which  is  df/dE. 

13 

E  is  the  energy  in  appropriate  units.  Berkowitz   gives  the 

relation  which  is 


a 

2 

ire  h 

mc 

df 
dE 

Comb: 

Lning 

Eqs. 

(1) 

to 

(5) 

yields 

n 

-  1  ■ 

2TTNe2    'S 

—^ 

f   1 
s 

r    2          2. 

[CO       -    CO    ) 

s 

(5) 


2    2,2    2  2 

(co  -  co  )  +  y  ^ 
s   s  s 


...    od    a(E2  -  E2)dE 
cNh  r    __js s 

„  2  /   ,2   r2,2  ,  Jl  2JI 
2tt  ^T  (E  -  E  )  +  tr  y  E 
Is  s 


(6) 


In  Eq.  (6),  c  is  the  velocity  of  light,  and  h  is  Planck's 
constant.   The  integral  has  been  written  in  terms  of  energy, 
E  ■  hco/2iT  =  tno.   Frequently  cross  section  data  are  given  in 

terms  of  energy  instead  of  frequency;  hence,  the  form  of 

» 
Eq.  (6)  is  more  convenient. 

When  y  is  equated  to  zero,  Eq.  (6)  reduces  to  Eq.  (1)  of 

14 
Liggett  and  Levinger.    This  check  provides  confidence  in  the 


algebraic  manipulation  involved  in  developing  Eq.  (6).   Retain- 
ing Y  in  Eq.  (6)  avoids  mathematical  difficulty  when  w  =  oo 
s  s 

or  E  =  E  . 
s 

Equation  (6)  is  correct  for  N  molecules  in  the  ground  state. 

If  the  gas  is  excited,  the  influence  of  excited  species  must 

be  included.  Excited  species  occur  in  excimer  lasers.  To 

develop  an  analogous  equation  to  Eq.  (6),  consider  the 

hypothetical  atom  shown  in  Fig.  1.   For  ease  of  writing,  y 

is  equated  to  zero.  Define  subscripts  as  follows:  i,  initial 

state;  f,  final  state;  u,  upper  level;  and  £,  lower  level. 

The  oscillator  strength  is  always  written  as  f.f»  i.e.,  the 

first  subscript  is  the  initial  state.  Consequently,  f  «  is 

emission,  and  f .   is  absorption.   In  accordance  with  convention 

and  the  need  for  compatible  equations,  f  .  is  negative,  whereas 

fp  is  positive.  See  Aller   or  Kemble.    Not  every  author 

uses  the  same  notation.   For  example,  Mitchell  and  Zemansky 

use  f,..  for  subscripts, 
fi 

Having  discussed  notation,  the  seven  terms  implied  by 
Fig.  1  are  written  explicitly 
2 


,        2iTe 
n  -  1  = 


m 


r    f01N0        .      f12Nl       +     f21M2  f10Sl 

2  22  22  22  2 

L^Ol  "  W         W12  "  W         a)21  -  W         w10  "  W 


chN.       r       a-dE.             chN        r        <j-.dE., 
4. 0     /  0     Oc         1    /  1     lc  .  ... 

+  ^rJ    ~i 2+~t~J    ~i 2   +  etc*        (7) 

2tt  I     E^     -  E^        2ttZ        I     Ef     -  EZ 

0      Oc  1     lc 


The  subscript  c  identifies  continuum;  c  occurs  in  the  integrals 
The  relation  between  emission  and  absorption  oscillator 
strengths  is  given  by  Aller   as 

u 
where  g.  is  the  degeneracy  of  the  lower  state.   Using  Eq.  (8), 
Eq.  (7)  can  be  rewritten  as 


(9) 


-  i  -  2v*2  y 

m          — 
1=0 

y            ik  i 

n  • 

Z-           2            2 
k=i+l      ik 

^        n     S 
s=0 

oo 

r      a  dE 
/          s     sc 

J         2            2 

I      EZ      -  EZ 
s     sc 

In  Eq.  (9),  i  is  a  lower  state,  and  k  is  an  upper  state. 
Obviously,  Eq.  (9)  involves  only  absorption  oscillator 
strengths.   The  quantity  I  is  illustrated  in  Fig.  1. 

The  similarity  between  the  double  sum  in  Eq.  (9)  and  the 
gain  equation  is  apparent.   The  similarity  will  now  be 
discussed. 
B.   Index  of  refraction  from  gain 

The  influence  of  the  lasing  species  upon  the  index  of  re- 
fraction can  be  determined  from  the  gain.  Define  u  =  2ojk/c  as 
the  gain  coefficient.   From  Eq.  (9),  the  index  of  refraction 
due  to  the  lasing  species  can  be  written  as 


n  -  1  - 


2          (to..    -  u))f  ..N. 
TTe ik  lk  x 

mto     ,  N2,,        /0\2 

(to       -  a))     +  (yik/2) 


"Nigk 


(10) 


To  obtain  Eq.    (10),    the  approximation 


2  2   ~ 

w.,    -  co     =  2oo(co.,    -  (u) 
ik  ik 


was  used.   Equation  (10)  is  for  the  lasing  transition, 
gain  coefficient  is 


(11) 

The 


y  - 


Tie 


'  ik  ik  i 


meg. 


1  - 


N,g. 
k  l 


(12) 


Tc  (ulk  -a0*+  (Y.k/2)     I        "i5k  J 
Forming  the  ratio  of  n  -  1   to  u  gives 

n  -  1        (wik  "  u)c8k       (a)ik  "  ^Xsk 


^ik 


2^ik 


(13) 


If  the  laser  is  tuned  to  the  line  center,  i.e.,  u).,  =  go,  the 
index  of  refraction  due  to  the  lasing  molecule  is  zero.  When 
co.,  does  not  equal  co,  Eq.  (13)  can  be  used  to  estimate  n  -  1 
given  u.   Equation  (13)  applies  to  the  case  of  natural  line 
broadening.   Other  equations  can  be  derived  for  other  line 
broadening . 

C.   Sensitivity  of  index  of  refraction  to  specie  concentration 
Replace  a  typical  term  in  the  summation  of  Eq.  (9)  with  the 
term  shown  in  Eq.  (10).  Multiply  the  term  in  square  brackets 
in  Eq.  (10)  by  N./g..   To  find  the  sensitivity  of  n  -  1  to 


species  concentration  N . ,  differentiate  yielding 


oo 


a_      2    (oj   -  u))g.f  ^   a  dE 

on  _  jre ik 1  ik      ch   r    s  sc      ,.. . 

s 
Equation  (14)  is  useful  in  determining  the  relative  importance 
of  a  particular  state  or  electronic  configuration,  i.e.,  the 
i-th  state,  of  a  specific  chemical  species. 
D.   Dependence  of  index  of  refraction  on  mass  fraction 

Based  on  the  discussion  to  this  point,  one  additional 
summation  should  occur  in  Eq.  (9).   As  written,  Eq.  (9)  is  for 
a  single  chemical  species.  The  index  of  refraction  should  be 
calculated  by  summing  over  the  contributions  of  each  chemical 
species.  The  third  summation  is  accomplished  by  replacing  N 
in  Eq.  (9)  by  (N./NJ  (N-/N)  .   The  quantity  N./N-,  which  is  the 
fraction  of  species  I   in  state  i,  reflects  the  influence  of 
populations.   N.  is  the  number  density  of  species  £;  and, 
hence,  N./N  is  the  mole  fraction  X-  of  species  £  in  the  gas 
mixture. 

As  written,  Eq.  (9)  provides  the  contribution  of  species  I 
to  the  index  of  refraction;  identify  that  contribution  by  the 
symbol  Tp .   The  index  of  refraction  for  a  gas  mixture  having 
X.  mole  fractions  and  populations  N./N.  is 

n  -  1  =  N  2XJlT£  (15) 

i 
The  total  number  density  is  N. 


To  ascertain  the  dependence  of  n  on  mass  density,  introduce 
the  partial  density  pp  which  is  N.m? ;  mp  is  the  mass  of  one 
molecule  of  species  I.      A  mass  fraction  Yp  =  Pn/P  can  be  in- 
cluded in  Eq.  (15).  When  the  changes  are  made,  Eq.  (15) 
becomes 

n  -  1 .  -  p  Y-^r-  (16) 

£  *- 

Equation  (16)  shows  that  n  -  1  is  proportional  to  p;  however, 

Eq.  (16)  further  demonstrates  that  n  -  1  depends  on  mass 

fraction  of  each  species.   Since  T„  involves  a  summation  over 

populations  [refer  to  Eq.  (9)],  n  -  1  also  depends  on 

populations. 

Based  on  the  preceding  discussion,  An/n  depends  not  only  on 
Ap/p  but  also  on  populations  and  mass  fractions.   A  better 
choice  for  specifying  the  constraint  on  optical  medium 
inhomogeneity  is  An/n. 

To  estimate  the  maximum  allowable  An/n,  assume  the  maximum 
allowable  phase  shift  at  the  output  aperture  is  A/20.   The 
geometrical  length  of  the  optical  path  in  the  laser  cavity  is 
L.  The  reference  index  of  refraction  is  n  .  Under  these 


conditions 


^  =  _Q =  — X (17) 

nQ     nQ     (20)  (L)  U/; 
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Inserting  typical  values  for  an  excimer  laser,  the  maximum 

—8 
allowable  An/n  is  2  x  10 

II.   BACKGROUND  DATA  FOR  CALCULATION  OF  INDEX  OF  REFRACTION 

A.   Typical  species  concentrations 

During  an  interval  extending  from  the  early  part  of  the 
pulse  to  the  end,  a  steady  state  condition  prevails  for  many 
of  the  chemical  species.   Lasing  occurs  during  the  steady 
state  period.   Table  I  gives  the  number  density  for  30 
chemical  species  for  an  e-beam  driven  excimer  laser  with 
initial  concentrations  Ne:Xe:NF  : :0. 992:0. 006:0. 002.   The 
results  shown  in  Table  I  were  calculated  using  an  excimer 
laser  kinetics  computer  program  developed  by  Johnson,  Palumbo, 
and  Hunter. 

The  chemical  species  in  Table  I  are  listed  in  order  of 
decreasing  number  density.   Initial  mass  density  of  gases  can 
be  inferred  from  the  number  densities.   The  information  shown 
is  for  a  time  of  250  nanoseconds  from  initiation  of  e-beam 
pulse.   At  250  nanoseconds,  steadv  state  has  been  attained. 

The  dominant  molecule  is  Ne  with  a  number  density  about  200 
times  greater  than  Xe.   The  index  of  refraction  for  the  gas 
mixture  should  be  dominated  by  Ne.  However,  other  chemical 
species  present  in  trace  amounts  may  contribute  significantly 
due  to  large  oscillator  strength  or  due  to  resonance  at  laser 
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TABLE  I.  Number  densities  for  thirty  chemical  species  in  a 
XeF  excimer  laser. 


Ne 

9.7E19 

Xe 

5.8E17 

NF3 

1.8E17 

NF2 

9.0E15 

F 

7.9E15 

F2 

3.0E15 

hv 

2.1E14 

f" 

9.7E13 

p* 

9.3E13 

Ne* 

8.3E13 

e 

4.9E13 

NF 

3.6E13 

Ne+ 

3 . 1E13 

NeXe+ 

3.0E13 

F* 

2.4E13 

Xe+ 

2.3E13 

Ne* 

2.2E13 

Ne+ 

2.0E13 

XeF* 

1.6E13 

XeF 

1.3E13 

NF+ 

9.9E12 

Xe* 

9.7E12 

< 

8.4E12 

NeF* 

5.3E12 

Ne* 

5.2E12 

Xe* 

1.7E11 

N 

9.3E10 

Xe* 

7.8E10 

NF+ 

1.7E09 

N  F 
W2r2 

4.3E06 

-3 

Note:  hv  is  the  photon  density  cm  "  for  wavelength  of  352  nm. 

Units  for  number  density  are  cm 
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wavelength.   Initial  efforts,  which  are  reported  in  this  paper, 

focus  on  Ne.   Fortunately,  considerable  spectroscopic  informa- 

18—35 
tion  exists  for  the  ground  and  excited  states  of  Neon. 

Since  Ne  is  the  lasing  molecule  in  the  ubiquitous  HeNe  laser, 

research  in  the  early  and  mid  1960's  was  directed  toward 

thorough  understanding  of  Ne. 

B.   Oscillator  strengths 

The  energy  level  diagram  for  Neon  is  shown  in  Fig.  2.   The 
energy  levels  are  shown  for  each  electronic  configuration. 
Figure  2  serves  as  a  reference  for  subsequent  tables  which  are 
arranged  in  the  form  of  matrices. 

Table  II  is  a  matrix  which  shows  the  references  used  to 
obtain  oscillator  strengths.   For  example,  data  for  the 
transition  3s  -*■  3p  can  be  found  in  references  18,  20,  22, 
26-29,  33,  and  34.   The  reason  that  abundant  data  exist  for 
3s  •+■  3p  is  that  the  HeNe  lasers  employ  this  transition.   When 
an  "X"  occurs  in  Table  II,  the  transition  does  not  occur  due 
to  selection  rules. 

Not  all  data  are  given  conveniently  as  oscillator  strength.' 
The  data  may  be  in  the  form  of  lifetime,  Einstein  coefficient 
A  0 ,  or  matrix  element.   All  data,  except  for  3s  ■+■  4p  transi- 
tion, were  converted  to  oscillator  strength  for  an  array. 
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Methods  for  changing  line  or  multiplet  f-numbers  to  array 

-1  C  O  f. 

values  are  given  in  Aller  "  and  Condon  and  Shortley. 

When  a  blank  occurs  in  Table  II,  the  oscillator  strength 

12  14 
was  estimated  by  extending  the  continuum      or  by  extrapola- 
tion of  plots  of  f-number  versus  principal  quantum  number,  n. 

The  3s  ■>  4p  transition  for  Neon  has  lines  in  resonance  with 
XeF  laser  radiation  in  351  to  355  nm  wavelength  interval. 
Consequently,  Eq.  (9)  must  be  evaluated  line-by-line.   Large 
values  of  n  -  1  may  occur  due  to  lines  in  the  3s  ■*■  4p  array. 

The  array  oscillator  strength  is  obtained  from  Gruzdev  and 

32 
Loginov,   while  the  relative  line  strengths  can  be  determined 

31 
from  Table  4  of  Semenov  and  Strugach. 

C.   Photoionization 

According  to  Eq.  (9),  photoionization  influences  the 

magnitude  of  n  -  1.  Data  are  available  for  the  photoioniza- 

14 
tion  cross  section  of  Neon  in  the  ground  state.    These  data 

were  used  to  evaluate  the  integral  in  Eq.  (9). 

2  2  5 
For  Neon  in  the  Is  2s  2p  3s  (metas table)  configuration,  the 

photoionization  cross  section  can  be  obtained  from  McCann  and 

37  38 

Flannery   or  from  Hazi  and  Rescigno. 

For  other  configurations  of  Neon,  information  about  C  could 

not  be  found.   One  approach  to  supply  the  missing  data  is  to 

calculate  a  using  procedures  outlined  in  references  13  and 

37-41.   The  calculations  for  O   were  beyond  the  scope  of  the 
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present  study.   The  contributions  of  a  due  to  higher  energy 

s 

configurations  of  Neon  have  not  been  included. 
III.   XeF  EXCIMER  LASER 

A.  Typical  conditions 

Four  aspects  of  excimer  lasers  are  relevant  to  beam  quality 
as  follows:   number  density  of  various  chemical  species, 
populations,  atomic  and  molecular  data  for  oscillator  strength 
and  photoionization  cross  sections,  and  the  flow  of  gases  in 

the  laser  channel.   Both  number  density  and  populations  in- 

17  42 
volve  a  detailed  knowledge  of  kinetics.   *    For  the  calcula- 
tions of  this  paper,  the  inhomogeneities  due  to  waves  are  not 
included.  Suppression  is  a  separate  and  severe  problem.  ' 
Table  I  provides  information  about  species.   Table  III 
gives  conditions  for  a  typical  XeF  laser  with  flow. 

B.  Transient  conditions 

Power  extraction  occurs  for  a  period  the  order  of  a  micro- 
second.  Variations  in  mass  density,  p,  occur  on  a  time  scale, 
T,  the  order  of 

T  -  f  (18) 

where  d  is  a  characteristic  length  and  a  is  the  acoustic 
velocity  which  is  300  to  400  m/sec.   Since  T  is  several 
hundred  microseconds,  changes  in  p  do  not  occur  during  the 
pulse  of  e-beam  and  lasing. 
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TABLE  III.   Typical  conditions  in  XeF  flowing  excimer  laser 

Initial  Gas  Composition,  Mole  Fractions 

Ne  95.5 

Xe  4.3 

NF3  0.2 

Initial  Gas  Temperature,   K  300 


Initial  Gas  Pressure,  Atmospheres  3-4 

Temperature  Rise  in  Cavity  due  to  Pulse,  K        100 

Pulse  Length,  Microseconds  0.5-2 

Flow  Mach  Number  0.5 

Another  factor  which  is  related  to  flow  is  turbulence. 

47 
Turbulence  can  affect  adversely  beam  quality.    The  influence 

of  turbulence  is  ignored.   Consequently,  the  results  of  this 

paper  are  applicable  to  a  laser  which  has  perfect  damping 

between  pulses. 

IV.   CALCULATION  OF  INDEX  OF  REFRACTION 

A.  Plasma  and  free  electrons 

The  contribution  to  the  index  of  refraction  of  free  elec- 

48 
trons  in  the  plasma   is 

a)2 
n2  =  1  -  -|  (19) 

where  oj  is  the  plasma  frequency  given  by 
P 
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2 

2   neS  4  7—  2 

uj  =  -2—  =  [5.63  x  10  VJT]  (20) 

P   £nm  e 

r    0  e 

-3 
In  the  third  term  of  Eq.  (20),  n  had  units  of  cm  .   As  an 

example,  from  Table  I,  n  =  4.9E13  giving  to  =  4E11.  At  a 

wavelength  of  352  nm,  the  value  of  w  is  5.4E15.  Hence  n  -  1 

due  to  electrons  has  a  value  of  3E-07.  When  n  approaches 

1E17,  the  optical  properties  of  the  free  electrons  become 

important. 

B.  Index  of  refraction  of  ground  state  Neon 

Using  Eq.  (6)  with  y     =0,  the  oscillator  strengths  from 
references  14,  21,  26,  32,  and  33,  the  photoionization  cross 

section  from  reference  14,  and  the  wavelengths  from  reference 

49 
24,  Etchechury  '  calculated  n  -  1  to  have  a  value  of  1.0006705 

at  352  nm.   The  value  compares  favorably  with  that  given  by 

Liggett  and  Levinger. 

Of  interest  is  the  size  of  the  photoionization  integral 

relative  to  the  sum  of  discrete  terms  in  Eq.  (6).   Table  IV 

summarizes  the  values  of  n  -  1  due  to  photoionization.   At 

352  nm,  the  photoionization  integral  contributes  about  9%  of 

the  overall  value  of  n  -  1. 

C.  Index  of  refraction  of  me tas table  Neon 

To  calculate  (n  -  1)  due  solely  to  Ne(3s) ,  Eq.  C.6)  is 
used.   The  summation  includes  both  emission,  e.g.,  3s  ■*  2p, 
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TABLE  IV.   Contribution  of  photoionization  to  n  -  1  for  ground 
state  Neon. 

Energy,  Rydbergs    Wavelength,  nm        (n  -  1)  x  10 

0.1  914  5.835 

0.2  457  5.864 

0.3  305  5.915 

0.4  228  5.988 

0.5  183  6.085 

0.6  152  6.211 

0.7  130  6.370 

0.8  114  6.568 

0.9  102  6.818 

1.0  91  7.134 


and  absorption,  e.g.,  3s  ■*   3p,  terms.   Only  one  integral  is 
necessary  for  photoionization  which  is  for  3s  ->  continuum. 
Equation  (9)  could  be  used  if  all  N.  or  N  except  N„  were 

zero. 

49 
Etchechury  '  calculated  n  -  1  due  only  to  Ne(3s) .   Table  V 

shows  the  data  for  calculating  (n  -  1)  due  to  discrete  transi- 
tions.  For  transitions  3s  -*■  2p,  3p,  5p,  the  individual  line 
transitions  are  replaced  by  an  equivalent  array.  However, 
for  the  transition  3s  +   4p,  each  line  must  be  accounted  for. 
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Some  lines  in  the  3s  -»■  4p  array  are  resonant  with  the  XeF 

photons. 

TABLE  V.   Information  for  discrete  transitions 


Transition 

Array  Oscillator 

Energy , 

Number  of  Terms 

Strength,  f 

Rydbergs 

3s  ■+  2p 

-0.0482 

1.230 

One  using  f . 

3s  •*■  3p 

1.080 

0.146 

One  using  f . 

3s  -*■  4p 

0.0106 

0.259 

Thirty  using  f 

3s  +  5p 

0.00129 

0.302 

One  using  f 

3s  ■*■  np 

Included  in 
continuum 
integral  for 
n  >  5 

Note:   f  means  oscillator  strength  for  an  array;  f  means 

A  JL 

oscillator  strength  for  a  line. 


For  the  transitions  3s  ■*  np(n  >  5)  ,  the  effect  is  included 
in  the  integral.   An  equivalent  cross  section  is  determined 
from  the  spectral  density  of  oscillator  strength.   References 
12  and  40  discuss  the  method. 

Table  V  gives  the  array  oscillator  strength,  f  ,  for 
3s  ■*-  4p  resonant  transition.  To  obtain  the  line  oscillator 
strength,  f  ,  for  each  of  the  30  lines,  the  table  of  relative 
line  strengths  from  reference  31  can  be  used.   The  results  are 
shown  in  Table  VI.   K  is  an  index  for  the  upper  levels  (4p) 
while  M  is  an  index  for  the  lower  levels  (3s) . 
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Using  the  information  from  Table  VI,  a  graph  of  (n  -  1)  due 
only  to  3s  -*■  4p  transition  was  calculated.   Figure  3  shows 
four  resonances  in  the  XeF  lasing  region.   Also  shown  in  Fig.  3 
are  the  wavelengths  for  the  line  centers  and  the  specific 
transitions.  For  the  specific  transitions,  Paschen  notation  is 
used. 

At  a  wavelength  of  352  nm  and  a  Ne(3s)  density  of 
2.68  x  10   cm  ,  Etchechury   calculated  (n  -  1)  due  solely 
to  Ne(3s) .   The  calculated  result  is  n  -  1  =  -  7.79E-10.   At 
the  same  wavelength  and  number  density  for  Ne(3s) ,  the 
contribution  of  photoionization  to  n  -  1  is  3.4E-12.   The 
percentage  contribution  of  photoionization  is  0.4%  at  352  nm. 

Later,  the  ratio  (n  -  1)/N  will  be  used;  the  ratio  is  the 
index  of  refraction  minus  unity  per  molecule.   In  terms  of 
this  ratio,  the  results  for  Ne(3s)  at  352  nm  are  as  follows: 
overall  n  -  1  is  -  2.907E-23  and  due  to  photoionization, 
+  1.27E-25. 

Table  VII  summarizes  the  contributions  to  (n  -  1)  of  the 
various  transitions.   The  table  has  been  calculated  for  a 
number  density  equal  to  Loschmidt's  number,  and  the  values  of 
the  wavelength  have  been  selected  for  XeF. 
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D.   Index  of  refraction  of  Neon  in  higher  energy  levels 

Using  Eq.  (6),  n  -  1  was  calculated.   Information  concern- 
ing the  photoionization  cross  section  for  energy  levels  above 
3s  could  not  be  found.   The  trend  in  a  between  ground  state 
Neon  and  metastable  Neon  is  a  large  reduction;  a  is  a  factor 
of  ten  less  for  Ne(3s)  than  for  Ne(2p  ) .   If  the  trend  in  the 
magnitude  of  a   continues,  the  contribution  of  the  continuum 
should  be  small. 

For  Ne(3p)  and  higher  energy  levels,  the  energy  of  XeF 
photons  occurs  above  the  ionization  threshold.  Hence,  the 
integration  of  the  continuum  has  both  a  positive  and  a 
negative  contribution,  i.e.,  the  denominator  of  the  integrand 
in  Eq.  (6)  switches  from  negative  to  positive  as  energy 
increases. 

In  Eq.  (6),  the  damping  term  was  equated  to  zero.   All 
transitions  between  excited  configurations  can  be  repre- 
sented in  a  matrix  as  shown  in  Table  VIII.   Each  element  in 
the  matrix  is  the  contribution  to  n  -  1  per  molecule  for  the 
transition  represented  by  the  matrix  element.   As  an  example, 


consider  3d  •*■   4p  transition  which  is  absorption;  assume 
3d 


-3 

N,  ,   is  3E11  cm   .   The  contribution  of  3d  -*■  4p  is 


(n  -  1)  =  (^-^N  x  10""27=  (-  7695)  (3  x  1011)  (10_27) 

-12 
=  -  2.3  x  10 
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Transitions  above  the  diagonal  in  Table  VIII  are  emission 
as  denoted  by  the  letter  E.   Transitions  below  the  diagonal  in 
Table  VIII  are  absorption  as  denoted  by  the  letter  A.  When  the 
columns  are  summed,  the  contributions  to  n  -  1  per  molecule 
for  all  transitions  are  included.   For  the  example  above,  the 
contribution  of  Ne(3d)  is 

(n  -  1)  =  (31324) (3  x  1011) (10~27)  =  +  9.4  x  10"12 
In  a  similar  manner,  all  transitions  can  be  calculated  and 
summed  if  the  populations  are  known. 

Table  VIII  was  calculated  for  a  wavelength  of  352  nm.   The 
results  are  not  sensitive  to  wavelength  except  for  the  two 
resonances  indicated  by  R*.   The  two  resonant  transitions 
are  between  3s  and  4p  either  in  emission  or  absorption. 

V.   CALCULATION  OF  NUMBER  DENSITIES  AND  POPULATIONS 

42 
Cole-   modified  the  computer  program  of  Johnson,  Palumbo, 

and  Hunter   to  yield  populations  of  excited  Neon.   The 

program  prior  to  Cole's  modification  calculated  Ne*,  which  is 

the  population  of  excited  Neon  irrespective  of  energy  level. 

Cole  divided  Ne*  into  two  groups:  metas table  Neon,  i.e., 

Ne(3s) ,  and  one  number  representing  all  higher  levels. 

As  an  example,  the  following  reaction  generates  Ne*, 

XeNe  +  e"  ■*■   Ne*  +  Xe  (21) 
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with  a  rate  coefficient  6E-08.   The  reaction  given  by  Eq.  (21) 
was  replaced  by  two  reactions 

XeNe+  +  e~  -*•  Ne(3s)  +  Xe  (22) 


and 


XeNe+  +  e~  +   Ne(z)  +  Xe  (23) 


In  Eq.  (23)  Ne(Z)  represents  Ne  in  any  energy  level  above  3s. 
The  reaction  rate  coefficients  for  Eqs.  (22)  and  (23)  were 

3.60E-08  and  2.40E-08  respectively. 

19 

The  volumes  by  McDaniel,  et  al.,    were  helpful  in  estimat- 
ing the  branching  ratios;  however,  many  of  the  choices  were 
arbitrary.   The  modified  computer  code  can  be  altered  to  re- 
flect improvements  in  rate  constants. 

To  distribute  the  Ne(Z)  populations  among  Ne(3p) ,  Ne(4p) , 
Ne(4s) ,  etc.,  a  Boltzmann  distribution  was  assumed.   The 
temperature,  which  varied  as  time  evolved,  was  determined  from 
the  ratio  of  populations  Ne(3s)/Ne(2p  ).   Knowing  the  tempera- 
ture, the  energy  level  diagram,  and  the  degeneracies  of  the 
levels,  the  Ne(Z)  could  be  distributed. 

When  the  populations  were  determined,  the  results  shown  in 
Table  VIII  were  used  to  calculate  (n  -  1)  for  each  electronic 
configuration.   Finally,  all  contributions  from  each  electronic 
configuration  were  summed  to  obtain  (n  -  1) . 
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The  calculations  were  performed  twice,  once  with  a  current 

2  2 

density  of  lOA/cm  and  again  with  llA/cm  for  the  e-beam.   The 

e-beam  voltage  was  constant.  The  energy  deposited  in  the  gas 

varied  by  10%. 

A  quantity  (n  -  l)/(n  -  l)n  was  determined  as  a  function  of 
time.   The  quantity  (n  -  l)n  is  the  value  at  the  start  of  the 
e-beam  pulse.   For  a  wavelength  away  from  any  of  the  3s  ■*■  4p 
resonances,  (n  -  l)/(n  -  l)n  changed  from  unity  at  zero  time 
to  0.999825  for  low  power  and  0.999810  for  high  power.  The 
fractional  change  is  1.5E-05.   In  terms  of  phase  shift,  6,  a 
laser  cavity  with  length,  L,  less  than  35  meters  will  have 
6 A  less  than  0.1. 

For  the  case  of  a  wavelength  near  a  resonance  in  the 
3s  -*■  4p  transition,  the  consequences  are  different.   A  length 
of  laser  cavity  of  only  0.1116  meters  yields  6/X  equal  to  0.1. 

Recall  the  power  variation  was  10%. 

42 
Details  of  the  calculations  are  given  in  Cole's  thesis. 

VI.   CONCLUSIONS 

Neon,  which  is  a  dominant  molecule  in  a  XeF  excimer  laser, 

has  transitions  in  the  3s  ■+■  4p  array  which  are  resonant  with 

laser  photons.   Operation  of  the  XeF  laser  near  one  of  the 

resonances  could  cause  extreme  sensitivity  to  small 

variations  in  power  density. 
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FIG.  1.  Energy  level  diagram  for  a  hypothetical  atom  used 
to  derive  equation  for  n  -  1. 
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FIG  2.      Energy  level  diagram  for  Neon.      (From  R.   J.   Pressley, 
Handbook  of  Lasers,   Chemical  Rubber  Co.,   Cleveland 
Ohio,   197170) 
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